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Some properties of the code newly developed to simulate the neutron 
detection process in a NMM are briefly described. 

PACS numbers: 29.40.Mc, 29.85. +c 



1. Introduction 

Recent ALADIN Collaboration experiments [1] presented an extensive 
evidence that the hot projectile prefragments created in the peripheral to 
semicentral relativistic heavy-ion collisions are thermally equilibrated. Tem- 
perature (T) and excitation energy (E x ) are the essential variables specify- 
ing their thermodynamic properties. The dependence of T on E x (caloric 
curve) revealed a plateau [2] interpreted as indication of the liquid-gas phase 
transition in finite nuclear systems. Therefore, a progress in measurement 
of E x and T is a necessary prerequisite for the further comprehensive study 
of the phenomenon. 

The predictions of Friedman [3] indicate that an excitation energy of 
a prefragment, whose mass and charge are generally unknown, can be de- 
termined by means of the multiplicities of the emitted neutrons and light 
charged particles (Z<2). We intend to apply this method to the target spec- 
tators by employing a large-volume Gd-loaded liquid scintillation detector 
(a neutron multiplicity meter - NMM) to count neutrons originating from 
their decay. A study of feasibility and limitations of the Friedman's method 
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in the ALADIN environment asked for a program giving a more compre- 
hensive description of the neutron detection process in a NMM than offered 
by DENIS [4], presently used in several laboratories [5]. These considera- 
tions prompted us to develop MSX, a Monte-Carlo code containing several 
improvements in comparison with its predecessor. 



2. Neutron detection process in a NMM 

The detection process in a NMM [5] consists of the following stages: 
i) Neutron slowing down due to the multiple scattering on the hydrogen and 
carbon nuclei of the scintillator. Ultimately, the neutron comes into thermal 
equilibrium with the medium and diffuses within the scintillator volume. 
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Fig. 1. Distribution of capture events within the volume of the tank for 1 MeV 
(left) and 30 MeV (right) neutrons. 

ii) Radiative capture following one of the three possibilities: in 1 H with the 
emission of a single 2.224 MeV 7-ray, in 155 Gd or 157 Gd accompanied by 
the emission of a cascade of 7-rays with the summed energies of 8.536 and 
7.937 MeV, respectively. 

Hi) Gamma-ray thermalization due to the consecutive Compton scattering 
from the electrons of H and C atoms of the scintillator. 

iv) Electron stopping in the medium. Part of the electron's energy loss, 
determined by the scintillation efficiency, is converted into visible light. 

v) Transport of light to the photocathodes of photomultipliers (PMs) through 
multiple reflections on the diffusely reflecting coating of the interior of the 
tank. 

The above stage (i) involves two vastly different time scales. The dissi- 
pation of neutron's total initial kinetic energy occurs within a nanosecond 
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time scale, giving the initial prompt part of the scintillation signal. The 
signal coming from the 7-rays follows with a random delay of up to 50 /is, 
the time scale determined by the probability of capture, terminating the 
diffusion process. Therefore essential for the operation of a NMM is the 
stretching of the capture probability over tens of microseconds which allows 
to count one-by-one over hundred neutrons emitted in the interaction. 



3. Novelties in MSX 

• Throughout the program 

A. Algorithm optimized to maximize the calculational speed. 

B. Memory utilization optimized with the aid of pointers allows to track 
simultaneously (practically) an arbitrary number of particles. 

• Stage (i) 

A. Use of the cross section files supplied by the international Nuclear Data 
Centers. They contain periodically updated data for many reactions on H, 

C. and Gd (23 total and 5 differential cross sections used by us). 

B. Exact treatment of aa/oYL for: ^(n^B, 12 C(n,n) 12 C, 12 C(n,n') 12 C 
(first, second and third excited state). This permits to follow the main 
tendencies of the angular distributions - rapid changes across the groups 
of prominent resonances in the n + 12 C interaction at low energies and in- 
creasing forward peaking as the energy exceeds the resonance region. 

C. Multiple neutron emission in the 12 C(n,2n) 11 C reaction taken into ac- 
count. It permits to determine the number of misclassified events (two 
neutrons detected instead of the incident one). 

• Stage (ii) 

A. Radiative capture in the resonance region for the 155 Gd(n,7) 156 Gd and 
157 Gd(n,7) 158 Gd reactions included. 

B. 7-ray cascades following neutron capture taken into account by combin- 
ing the available experimental data with the calculated continuous part of 
the spectrum. Calculations were done on base of the Hauser-Feshbach sta- 
tistical model for the decay of the J w = 2~ capturing state. Besides mean 
multiplicities, defining mean capture 7-ray energies, also fluctuations around 
the mean taken into account. 

• Stage (v) 

A. Simulation of light transport to the PM photocathodes incorporated into 
the program. It permits to determine the effect of electronic thresholds on 
the detection efficiency. Moreover, the distribution of photon arrival times 
determines the PM current pulse shape, thereby defining the double-pulse 
resolution. 

B. Calorimetric mode (in the stage of implementation). It will permit to 
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simulate the dependence of the prompt peak on the energy deposited in the 
scintillator by the incident neutrons. 

4. Representative results for the selected geometry 

The geometry of a NMM destinated to count neutrons from the target 
residue in coincidence with the charged products of the projectile-target 
interaction is shown in Fig. 1. The escape cone with 50° opening angle 
enables undisturbed detection of the intermediate mass fragments from the 
projectile-like source by means of the ALADIN heavy fragment detectors. 
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Fig. 2. Energy dependence of the neutron detection efficiency in the energy range 
0.5-32 MeV. 

Light charged particles arising from the intermediate-velocity source are 
registered with the 92-clcmcnt Si-CsI(Tl) array [1]. These are used to tag 
the impact parameter involved in the collision. In the same time neutrons 
from the latter two sources are prevented from interacting with the scintil- 
lator and the tank walls. The results of MSX simulations for the indicated 
geometry are presented in Figs. 1 and 2. 

REFERENCES 

[1] A. Schuttauf et al, Nucl. Phys. A607, 457 (1996). 
[2] J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995). 
[3] W.A. Friedman, Phys. Lett. B242, 309 (1990). 

[4] J. Poitou, H. Nifenecker, C. Signarbieux, Report CEA-N-1282 (1970); 

J. Poitou, C. Signarbieux, Nucl. Instr. and Meth. 144, 113 (1974). 
[5] J. Galin, U. Jahnke, J. Phys. G: Nucl. Part. Phys. 20, 1105 (1994). 



